HD-THEP 99-41 
September 1999 



THE CP- VIOLATING TRIPLE GLUON 
INTERACTION IN Z -> 4 JETS 1 

O. Nachtmann and C. Schwanenberger 

Institut fur Theoretische Physik, Universitat Heidelberg, Philosophenweg 16 

D-69120 Heidelberg, Germany 
O.Nachtmann@thphys.uni-heidelberg.de, C.Schwanenberger@thphys.uni-heidelberg.de 



ABSTRACT 

We analyse CP-violating effects in Z — > 4 jet decays, assuming the presence 
of a CP-violating effective triple gluon coupling. We discuss the influence of 
this coupling on the decay width. Furthermore, we analyse different CP-odd 
observables and propose strategies of a direct search for such a CP-violating 
GGG coupling. The present data of LEP 1 should give significant information 
on the coupling. 



1 Supported by German Bundesministerium fur Bildung und Forschung (BMBF), 
Contract Nr. 05 7HD 91 P(0), and by the Landesgraduiertenforderung 



1 Introduction 



In electron-positron collider experiments at LEP and SLC, a large number of Z bosons has been 
collected so that the detailed study of the decays of the Z boson has been made possible [j]. 
An interesting topic is the test of CP symmetry in such Z decays. There is already a number 



of theoretical (@-[19| and references therein) and experimental [20|-[27| studies of this subject. 
In the present paper we will study a flavour-diagonal Z decay where CP-violating effects within 
the Standard Model (SM) are estimated to be very small Thus, looking for CP violation 
in such Z decays means looking for new physics beyond the SM. 

For a model-independent systematic analysis of CP violation in Z decays we use the effective 
Lagrangian approach as described in Here we consider the 4 jet decays of the Z boson. In 
TRITBipSI the effects of CP violating couplings involving heavy quarks were studied for 3 and 4 



jet decays of the Z boson. On the other hand the 4 jet decays offer also the possibility to study 
the CP violating triple gluon coupling which was listed in p5| , |3T| and was investigated in pT| 
The following three subprocesses contribute to the 4 jet decay: 

: + (p+,v) e~ (p-,w) -» Z(p,j) -> q'(k-,s,B) q'(k + ,r,A) G{k 1 ,K 1 ,a) G(k 2 ,K 2 ,b) , 



e + (p+,v) e (p-,w) -> Z(p,j) -> q'(k_,s,B) q 1 (k +} r,A) q' (q_,u,D) q'(q + ,t,C) , 



(2) 



e+ (p+,v) e (p-,w) Z(p,j) -> q'(k_,s,B) q'(k +} r,A) q(q_,u,D) q(q +} t,C) , 

(Q*<f), (3) 

where q and q' denote quarks with flavour q, q' = u, d, s, c, b. We will always assume unpolarized 
e + , e~ beams. Only for process (fj) the CP- violating GGG coupling comes into play. We will 
show the results for this process alone as well as the results for the sum of them. In the 
experiments, of course, only the sum of the three processes can be observed easily. 

In chapter 2 we explain the theoretical framework of our computations. Next, in chapter 3, 
we analyse the anomalous coupling for partons in the final state. First, we discuss anomalous 
contributions to the decay width. Then, we investigate different CP-odd tensor observables as 
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in and calculate their sensitivities to the anomalous coupling. In order to find out how 
"good" for the measurement of the new coupling our observables are, we compare them to the 
optimal observable. In chapter 4 we study the optimal observable in a realistic scenario for an 
experimental analysis. Our conclusions can be found in chapter 5. 



2 Effective Lagrangian Approach 

For a model independent study of CP violation in 4 jet decays of the Z boson we use the 
effective Lagrangian approach as explained in ||. We add to the SM Lagrangian C SM the 
following CP-violating term which contains a mass dimension d < 6 local operator involving 3 
gluons: 



C CP (x) = ikS P { G aP {x) G^{x) G pc {x) } e a ^ g™ , (4) 

where G a v {x) represents the field strength tensor of the gluon and e a ^ fip is the totally anti- 
symmetric tensor with eoi23 = +1- A typical process where we find the corresponding vertex 
following from Cqp is shown in figure |l|. 




q 



Figure 1: Diagram for Z — > qqGG with the CP-violating vertex ([|). 
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We define a dimensionless coupling constant k using the Z mass as the scale parameter by 



k . (5) 



Here g s is the gauge coupling constant of QCD. For numerical calculations we set mz = 
91.187 GeV and a s = g 2 /4:7c = 0.118 since we consider a process at the Z mass scale [33 



Our calculations are carried out in leading order of the CP-violating coupling of Cqp and the 
SM couplings. All quark masses are neglected.^ 

3 Study of the CP-violating coupling for partons in the 
final state 

In this chapter we discuss an ideal experiment where one is able to flavour-tag the partons and 
measure their momenta. We present a study of the CP- violating coupling for process ([[]) with 
q' = u,c and q' = d,s,b separately and for the sum of the processes (0) ~ ©• We have computed 
the differential and integrated decay rates using FORM m and M Q for the analytic and 



VEGAS [42|] for the numerical calculation. We write the squared matrix element for the process 



(jl]) with final state p, 

p = uuGG, ddGG, ssGG, ccGG, bbGG , (6) 

in the form: 

R^iP) = 5 O (0)(P) + kSM {P) + k 2 S 2 (4>)^ . (7) 

Here stands collectively for the phase space variables, So denotes the SM part. For the 
processes (0), (Q) with final state p, 

p = uuuu, dddd, ssss, cccc, bbbb, 
uucc, 

uudd, uuss, uubb, ccdd, cess, cebb, 

ddss, ddbb, ssbb , (8) 
the matrix element contains only a SM part: 

R{4>)^ = Soi^ . (9) 



*For further details of the calculation we refer to 13£ 
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The definition of a 4 jet sample requires the introduction of resolution cuts. We use JADE 
cuts E31 requiring 



= 2^(1^) > 
mf 

with the angle between the momentum directions of any two partons (i ^ j) and Ei, Ej 
their energies in the Z rest system. The expectation value of an observable 0(<p) is then defined 

as 

fR(4>)d(!> ' 1 j 

3.1 Anomalous contributions to the decay widths 

The solid curves in figure |2] show the results of our calculations for the SM decay widths T SM 
as function of the jet resolution parameter y cut for process flip with the different final states p 
of (El). 




Figure 2: The decay width for different subprocesses as function of the jet resolution 
parameter y cut ([IZj/ The results for Z — ► ccGG (ssGG, bbGG) are identical to those 
ofZ^uuGG (ddGG). 
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To check our calculations we computed T also with the program COMPHEP and 
found — within numerical errors — complete agreement. 



Z -> 4 jets 




l 
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Figure 3: The 4 jet decay width (sum of processes - (j^) ) as function of the jet 
resolution parameter y cut (\T(\). 



As the decay width is a CP-even observable the contribution of the CP- violating interaction 
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to it adds incoherently to the SM one: 

r = r 5M + Ar CP , (12) 

with Ar CP being quadratic in the new coupling. In figure ^| the dashed curves represent AT CP 
as function of y cut assuming k — 1. 

In Figure || we compare for the sum of the processes ([!]) - (Q) T SM and Ar CP assuming 
K — 1. AT CP is only a correction of less than a per mille to r . Thus, considering the 
theoretical uncertainties in the SM 4 jet decay rate, a determination of the new coupling by 
measuring the decay width alone does not look promising. 

3.2 CP-odd observables for Z bbGG and Z ccGG 

Tagging of jets according to quark flavour or gluon is a difficult experimental task. What can 
be done routinely now is tagging of heavy quark b and c jets. Thus we study in this section the 
information obtainable on k (|j) from reaction ([!]) with q' = c and b. 

Let us first consider CP-odd observables constructed from the momentum directions of the 



q and q quarks (q = c, 6), k 9 = k 9 /|k 9 | and k ? = kg/|k^| (cf. flj^, pl]jr^j28|1 ) : 

= (k q - k q ) t (k 9 - x k^ |k, x k q \ n - 2 + (i ~ j) , (13) 

with i, j the Cartesian vector indices in the Z rest system and n = 1, 2, 3. 

The observables T^ 1 transform as tensors. For unpolarized e + e~ beams and our rotationally 
invariant cuts ([TOD their expectation values are then proportional to the Z tensor polarization 
Sij. Defining the positive z-axis in the e + beam direction, we have 



(S {j ) = - -10 . (14) 




This shows that the components are the most sensitive ones. 

Note that the tensor observables do not change their sign upon charge misidentification 
(k 5 <-> kg). Thus there is no need of charge identification in a measurement. We have also 



investigated vector observables like in pq] , but found them to be scarcely sensitive on the 
CP- violating coupling (Hf). 
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We have computed the expectation values of the observables (|T3|), for different JADE cuts 
as function of k. The expectation value of a CP-odd observable O has the following 
general form: 

YSM 

<0>= ck -jr- , (15) 

where T SM and T denote the corresponding Z — > 4 jets decay widths in the SM and in the theory 
with SM plus CP-violating coupling, respectively. In an experimental analysis T SM should be 
taken from the theoretical calculation, T and <0> from the experimental measurement. The 
quantity < O > ■ T is then an observable strictly linear in the anomalous coupling. 

The relative statistical error 5k to leading order in the anomalous coupling in a measurement 
of the coupling k using the observable O is given by: 



6K = ^T^r ' (16) 



where N is the number of events within cuts. A measure for the sensitivity of O to k is 1/ 5k. 

In addition to the tensor observables (|T3|) we study the optimal observable, which has 
the largest possible statistical signal-to-noise ratio [45j46|,47l. Neglecting higher orders in the 



anomalous coupling the optimal observable for measuring k is obtained from the differential 
cross sections (0) and @, respectively, as 

O(0) = f§ . (17) 
The expectation value has then the following form: 

<0>=ck, (18) 

with the coefficient 

S ™-S x {(j))d(i) . (19) 



/ 5*0^ J S Q {4>) 
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Figure 4: The inverse sensitivities of tensor and optimal observables to k obtainable in 
the subprocesses ^) with q' = c and q' = b as function of the jet resolution parameter 
Ucut ftTty ) assuming (fflty for the number of events. 
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We have calculated the sensitivities to k for different tensor observables and the optimal 
observable varying the jet resolution parameter y cut . We assume a total number of 5 • 10 5 4 jet 
events from ([!] - |3|) for y cut = 0.01: 

N(y cut = 0.01) = 500000 . (20) 

The number of events for other values of y cut and for the various subprocesses is then calculated 
within the SM. The total number of Z decays corresponding to ( p0[) is N to t — 6 • 10 6 . 

In figure | we show the inverse sensitivities to the CP-odd triple gluon coupling as calculated 
from (|T6|). The differences due to the different weight factors for tensor observables T33 (n = 
1,2,3) are only small but all these observables have significantly lower sensitivities than the 
optimal observable. Information on the spin of the final state partons in (|1]) - (|3]) is hardly 
available experimentally. Thus, we consider as observables only the parton's energies and 
momenta. Then, we didn't find any simple CP-odd observables with a significantly better 
sensitivity than those of the tensor observables T^ 1 (n = 1, 2, 3). It is thus of advantage to use 
the more complicated but much more sensitive optimal observable for the experimental analysis 
of the CP-odd coupling. 

In tables p] [| in appendix || we list the coefficient of the expectation value ([Tj|) for the 
optimal observable fli~7D for different values of the jet resolution parameter y cut (|i~(f) for the sum 
of the reactions (|l|) - (||) and for the reaction (H) with q' = c and b, respectively. 

4 CP- violating observables for untagged jets 

In this section we consider an experimental analysis of untagged jets which are ordered according 
to the magnitude of their momenta q«, i — 1, 2, 3, 4 : Q 

|qi| > |qal > > hil ■ ( 21 ) 

Thus, all processes (|TJ) - (^) are lumped together here. 

5 This we called analysis 4 in [ p8[ . Further details can be found therein. 
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Figure 5: The error (inverse sensitivity) <5k obtainable from the optimal observable 
as function of the jet resolution parameter y cut flip] ) assuming for the number of 



events. The solid curve represents the results for momentum ordered jets The 
dashed curve gives the sum of the results from reaction (j7[j with q' = c and b (cf. 
section \3.2J . The dotted curve represents the sum of the results from the reactions 
- (||) under the assumption that one is able to flavour-tag all partons and measure 
their momenta. 



The contributions to the decay width are as in the parton case in figure |[ Here suitable 
CP-odd tensor observables are 

T'\f = (qi - q 2 ), (qi X q^. |q x x q^"" 2 + (i » j) , (22) 

where n = 1, 2, 3 and qz = qz/|qz|- As expected from the results of section ||| these are by far less 
sensitive than the respective optimal observable for momentum ordered jets. For this reason, 
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we investigate in the following only the optimal observable for momentum ordered jets. 

We computed the inverse sensitivities 5k for the optimal observable ( |TTD for momentum 
ordered jets ( |2"T| ) using (|T^) as a function of the jet resolution parameter y cut (|TUp. The results 
are shown in figure || (solid curve). The loss of information about the parton charge and flavour 
due to the ordering of the jets according the magnitude of their momenta leads to a loss of 
sensitivity of the optimal observable. The difference between the dashed curve (results from 
reaction ([!]) with q' = c and b) and the dotted curve (results from the reactions ([!]) - @ 
under the assumption that one is able to flavour-tag all partons and measure their momenta) 
is essentially due to the different number of events as one can see from tables 0-0. 

In table | in appendix ^| we list the coefficient of the expectation value (|18|) for the opti- 
mal observable (|I7]) for momentum ordered jets (pl|) for different values of the jet resolution 
parameter y cut ([10]). 

5 Conclusions 

We have studied how one can search for CP violation in the 4 jet decays of the Z boson assuming 
a CP-violating triple gluon coupling, which can arise at one loop level in Higgs extensions of 
the SM |35| , in supersymmetric models |33 ] or in left-right symmetric models . 



We found that, for reasonable values of the coupling constants, the additional contribution 
of the contact interaction to the decay width is at most at the per mille level. The decay width 
alone is therefore not appropriate for determining the coupling constant. 

We have investigated different tensor observables as well as the optimal observables which 
can be used for the measurement of the anomalous coupling. The tensor observables have only 
a low sensitivity on the CP-violating coupling. No simple observables with sensitivities near to 
the optimal have been found. 

If it is possible to tag the flavour and the charge of all jets then, with a total number of 
Ntot(y C ut = 0.01) = 6 • 10 6 Z decays and choosing a jet resolution parameter y cut = 0.02 the 
anomalous coupling can be determined with an accuracy of order 0.3 at 1 s. d. level using the 
optimal observable. If flavour tagging is possible for b and c quarks and taking together the 
subprocesses (||) with q' = c and q' = b an accuracy of order 0.5 at 1 s. d. level is reachable. If 
no flavour is tagged and instead all four jets are ordered according to the magnitude of their 
momenta, the coupling constant k can be measured with an accuracy of order 1.5 using the 
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same total number of Z decays (see figure ||). 

In our theoretical investigations we assumed always 100% efficiencies and considered the 
statistical errors only. But the total number of Z decays collected by the LEP and SLC experi- 
ments together is of order 10 7 . Thus the accuracies in the determinations of k discussed above 
should indeed be within experimental reach. 



As shown in [32| the CP-violating triple gluon coupling can also be studied in the reaction 
pp — > 3 jets + X, for instance at the Tevatron. The accuracies obtainable there for k are 
generally of similar order of magnitude as for Z decays (cf. eqs. (4.10) and (4.12) of ||32|| ). In 
detail they depend, of course, on the number of available events within cuts. 

It is interesting to compare these accuracies to the limit on |ac| one receives under the 
assumption that the triple gluon operator (|j) delivers the main contribution to the electric 
dipole moment (EDM) d n of the neutron. At the moment the upper limit on d n is at 90% c. 1. 

m 

d n < 0.97- 10" 25 ecm . (23) 



Using the "naive dimensional analysis" described in |JT], but using the correct anomalous 
dimension for the 3-gluon-operator (|) as in f37j, we get 



\k\ < 2.2 • 10~ 7 . (24) 

This is much smaller than the accuracies reachable in our analyses with presently available 
numbers of Z decays. However, one should keep in mind, that many CP-odd operators can 
contribute to d n and cancellations among them cannot be excluded. 

To summarize: We have shown that the measurement of the optimal observables for flavour 
tagged and momentum ordered jets in 4 jet decays of the Z will give useful limits on k. (A 
FORTRAN-program for the optimal observables is available from the authors.^) Even if the 
limits obtainable are much worse than from the neutron's EDM they have the advantage of 
being direct, i. e. of involving explicitly 3 gluons in the splitting of one gluon into 2 gluon jets. 
Such studies with the LEP1 data could also serve as pilot studies for future investigations of 
this sort in the process e + e~ — > 4 jets at planned linear colliders (see e. g. |49[|), where due to 
the higher c. m. energy of 0.5 to 1 TeV effects of the dimension 6 effective 3 gluon coupling @ 
are of course enhanced. 



6 World Wide Web address: http://www.thphys.uni-heidelberg.de/~schwanen 
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Appendix A Numerical Values 

We list some numerical results for the coefficient of the expectation value for the optimal 
observable for partons and momentum ordered jets in the final state. The statistical errors of 
the numerical calculation are typically at the per cent level. 



Vcut 


c 


0.01 


3.014- 10- 5 


0.02 


6.158- HT 5 


0.05 


1.587- 10" 4 



Table 1: The numerical values of the coefficient of the expectation value ( 17^ for the 
optimal observable O = O ftTT{ ) for different values of the jet resolution parameter y cut 
(jTTJ) for partons in the final state (sum of the reactions ([][) - (TJjj. 



Vcut 


c 


0.01 


3.172 • 10~ 5 


0.02 


6.489 ■ 10~ 5 


0.05 


1.693- HT 4 



Table 2: The numerical values of the coefficient of the expectation value (fJ^j for the 
optimal observable O = O fiT\) for different values of the jet resolution parameter y cut 
f \n\ ) from reaction Z — > ccGG ([/[). 
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Vent 


c 


0.01 


3.159- lO" 5 


0.02 


6.515- 10~ 5 


0.05 


1.695- 10" 4 



Table 3: The numerical values of the coefficient of the expectation value ^TEj ) for the 
optimal observable O = O flH[) for different values of the jet resolution parameter y cut 
( jTZp from reaction Z — > bbGG 



Vcut 


c 


0.01 


1.709 • HT 6 


0.02 


2.666 • 10~ s 


0.05 


3.207- HT 6 



Table 4: The numerical values of the coefficient of the expectation value ftlBi ) for the 
optimal observable O = O ftP7{ ) for different values of the jet resolution parameter y cut 
( |77^ ) for momentum ordered jets (chapter^). 
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